Muon capture by ^He nuclei followed by proton and deuteron production 
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The paper describes an experiment aimed at studying muon capture by ^He nuclei in pure ^He 
and D2 + ^He mixtures at various densities. Energy distributions of protons and deuterons produced 
via fi~ + ^}ie ^ p + n + n + Ufi and fi~ +"^116 d + n + u^ are measured for the energy intervals 10 — 
49 MeV and 13 — 31 MeV, respectively. Muon capture rates, X^^p{AEp) and X'ii^p{AEd) are obtained 
using two different analysis methods. The least-squares methods gives X^^^ — (36.7 ± 1.2) s~^, 
Alp = (21.3 ± 1.6) s-\ The Bayes theorem gives A^p = (36.8 ± 0.8) s-\ AJ*^p = (21.9 ± 0.6) s'K 
The experimental differential capture rates, dXcs,p{Ep) /dEp and dXcs_p{Ed) / dEd, are compared with 
theoretical calculations performed using the plane-wave impulse approximation (PWIA) with the 
realistic NN interaction Bonn B potential. Extrapolation to the full energy range yields total proton 
and deuteron capture rates in good agreement with former results. 

PACS numbers: 34.50.-s, 36.10.Dr, 39.10.+j, 61.18.Bn 
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I. INTRODUCTION 

The study of few-nucleon systems is interesting and 
very important. It gives a microscopic description of 
complex systems within the framework of modern con- 
cepts of nucleon-nucleon interaction Using the nu- 
clear muon capture to study few-nucleon systems is a 
perfect tool since the nuclear structure had been found 
to play an important role in such systems 0, En- 
ergy transferred to a nucleus when muon capture oc- 
curred causes the excitation of low-lying levels in the 
residual nucleus up to the giant resonance region 0| or 
emission of intermediate-energy neutrons |^. This pic- 
ture is clear within the framework of the plane-wave im- 
pulse approximation (PWIA) J3_(and references therein) . 
However, some experiments Ul3, H, indicate that the 
energy transferred to the residual nucleus in muon cap- 
ture is large. It was found in those experiments that 
collective nuclear excitations like giant resonances play 
a decisive role in the muon capture process. In most 
cases the decay of the giant resonance was followed by 
the emission of a neutron and the formation of a daugh- 
ter nucleus in the above-threshold state for which it was 
then "beneficial" to decay via the proton or deuteron 
channel 0, S H IH 113 . 

An interesting feature of such nuclear decays is the 
emission of high-energy (40 — 70 MeV) charged particles 
(protons, deuterons) [ii, lH, QIIE EE 113 ■ By studying 
such an emission resulting from nuclear muon capture it 
is possible to get information both on the nuclear struc- 
ture and the muon capture mechanism itself 0, Q ■ The 
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emission of high-energy protons and deuterons in muon 
capture seems to be due to the existence of initial or 
final-state nucleon pair correlations and to a contribu- 
tion to the interaction from the meson exchange currents 
(MEG) [Hill. Note that the MEG contribution is very 
sensitive to the details of the wave function for the nu- 
clear system. 

In the region of large energy transfer (extreme kine- 
matics case) the MEG contribution to the interaction be- 
comes substantial. Note that MEG and nucleon-nucleon 
correlation effects are included "automatically" . For ex- 
ample, the calculation of the rate for muon capture by 
a deuteron poL |^ indicates that inclusion of MEG in 
the muon capture matrix element considerably increases 
the calculated capture rate at the boundary of the kine- 
matic region as compared to the contribution from the 
high-momentum components of the deuteron wave func- 
tion. The above mentioned factors may cause nuclear 
transitions with a large energy transfer. 

Though yields of charged particles in the muon cap- 
ture process are relatively small, the study of these events 
may give more information than other methods: it pro- 
vides an insight into the mechanism for excitation and 
decay of nuclei upon muon capture. So far, there is no 
microscopic description of the nuclear muon capture pro- 
cess . To ensure a correct comparison between theory 
and experiment, it is necessary to study muon capture in 
few-nucleon systems {A < 3), where a microscopic cal- 
culation of wave functions in the initial and final states 
is possible 20, 21]. 

Matrix elements calculations for the nuclear muon cap- 
ture transitions are usually performed using the wave 
functions model of the initial and final states. The wave 
function parameter values are chosen such that calculated 
and experimental data agree correctly for the case of 
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low-lying nuclear states spectra and corresponding mag- 
netic moments 0- In the case of light nuclei a multi 
particle shell model is frequently used. This model de- 
scribes (with a defined accuracy) these characteristics, 
i.e., the spectra and magnetic moments. However, the 
shell model accuracy may become insufficient because of 
poor knowledge of muon-nucleon interaction constants. 
In addition, there remains the problem of MEC. 

At present, general properties of nuclear transitions to 
the continuous spectrum for muon capture are treated on 
the basis of a resonant collective mechanism for the muon 
absorption by a nucleus 0, The strongest El tran- 
sitions, much like nuclear photo-disintegration reactions, 
form a giant dipole resonance and are collectivized into a 
continuous spectrum at muon capture '5| . The character 
of collective motions excited in nuclei at muon absorption 
is different from that in nuclear photo-disintegration re- 
actions. 

The giant resonance at muon capture differs from the 
photo-nuclear giant resonance by a greater importance 
of spin waves (similar to collective excitations in solids) 
and by a larger momentum transferred to the nucleus 
(neutrino momentum) for muon capture than for photon 
absorption with an energy in the vicinity of the giant res- 
onance. In addition, high-multipolarity transitions play 
a more significant part in muon capture than in photo- 
nuclear reactions. It is not yet clear why the charged 
particle yield at muon capture increases as one goes from 
Ip— shell nuclei to (2s — Id)— shell nuclei. Structure pe- 
culiarities of th e g iant resonance in (2s — Id)— shell nu- 
clei [23, lU 112, may play an important role, though. 

For example, the entrance states of one particle-one 
hole {Ip — Ih) nuclei should quickly decay into more 
complicated configurations which may emit various par- 
ticles before a thermodynamic equilibrium is established 
in the nucleus. This is the so-called decay from the pre- 
equilibrium state |3. l23|. In accordance with it, energy 
spectra of emitted protons and deuterons from {2p — 2h) 
states of the daughter nucleus must be well extended into 
the high energy region. In Ref. 0| the authors assumed 
that proton emission at muon capture may indicate the 
presence of {2p — 2d) states in the giant dipole configu- 
ration. 

While in the low-energy region of emitted charged par- 
ticles the resonant muon capture mechanism dominates, 
in the high-energy region the direct muon capture by 
correlated nucleon pairs seems to become prevailing. In 
the light of the aforesaid it is interesting to study muon 
capture by "^He (and '*Hc) nuclei followed by emission of 
protons 

/i^+^Hc^p + n + n + z/^, (1) 

and deuterons 

/i" +^Hc^d + n + i/^. (2) 

Note that muon capture by '^He is predominantly (70% 
of the cases) followed by the emission of tritons 

^~ -|-''^He -^t + Uf,. (3) 



However, this reaction was not studied in our experiment. 

Reactions and Q also attract interest because they 
are background reactions for the nuclear fusion process 
in the d/i'^He molecule 

d/^^He -^p + a + (4) 

to which considerable experiment alpl I2II27L I2ll29ll30| 
and theoretical |^ Is^, IH, 0, l35l| studies have been 
devoted in the past five years. In addition, the study 
of such systems gives the possibility of verifying funda- 
mental symmetries in strong interactions (such as charge 
symmetry or isotope invariance) js^H^ and solving some 
astrophysical problems . 

There has been only one experiment |39l |43 | in which 
differential probabilities for muon capture by ^He nuclei 
with the production of protons dX^^^/ dEp and deuterons 
dX^^p/dEd were measured at a few proton energies Ep 
in the range 17 — 52 MeV and deuteron energies Ed in 
the range 20 — 28 MeV. In addition, total summed rates 
for processes shown in Ens. Q and ^ were measured in 
three experiments j^^MEJand calculated in Refs. 0, 
2^23 . A recent review 0] is devoted to the experimental 
and theoretical study of the nuclear muon capture and in 
particular to the muon capture by He nuclei. It contains 
essentially the full list of theoretical and experimental 
work performed in this field until today. 

Other points indicating the importance and the neces- 
sity of studying processes of muon capture by '^He nuclei 
are the following: 

- Progress in the wave function calculations for the 
initial and final states of such a three-body sys- 
tem 47, 48, 49, 50, 51, 52] will give a better com- 
parison between experiment and theory. 

- Precise information on the characteristics of reac- 
tions ^ and Q in a "softer" proton and deuteron 
energy region as that in m 13 by using different 
techniques will be obtained. 

The purpose of the study described in this paper 
is to measure the energy distributions of protons and 
deuterons {S{Ep), S{Eii)) produced in reactions Q 
and ^ . We will also study the energy dependence of the 
differential probabilities for muon capture by '^He nuclei. 

II. EXPERIMENT 
A. Experimental set— up 

The experiment was carried out at the /iE4 chan- 
nel at the Paul Scherrer Institute (PSI) in Switzerland. 
The apparatus was originally designed and used to mea- 
sure the nuclear fusion rate in the molecular system 
dfM^Ee m mm 113. Figure □ schematically displays 
the apparatus as seen by an incoming muon. 

The cryogenic gas target, described in detail in [29l |. 
consisted of a vacuum isolation region ("V" in Fig.^ and 
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10 cm 



FIG. 1: Apparatus used in the /iE4 area. The view is that 
of the incoming muon. Note that the Tl and TO scintiUators 
are not shown. The labels are explained in the text. 



a cooled pressure vessel made of pure aluminum ( "T" in 
Fig. ^ . The pressure vessel enclosed a 66 mm diameter 
space which was filled with either pure ^He or D2 + '^He 
mixtures. Five stainless steel flanges held kapton win- 
dows over ports in the pressure vessel to allow the muons 
to enter and the particles of interest to escape from the 
central reaction region. In total, the target gas volume 
was « 250 cm'^. 

The incident muons, 17 x 10"^ /i/s at momenta 
34 MeV/c or 38 MeV/c were detected by a 0.5 mm thick 
plastic scintillator of area 45 x 45 mm^, called Tl, lo- 
cated at the entrance of the chamber. The electron im- 
purities in the muon beam were suppressed by a detector 
and a lead moderator, called TO, both having aligned 
= 44 mm holes, slightly smaller than Tl. Detectors 
TO and Tl are not shown in Fig. ^ since they lie in the 
plane of the paper. To reduce background coming from 
muons stopping in the entrance flange with their subse- 
quent nuclear capture and production of charged compo- 
nents (protons, deuterons, etc.), a 1 mm thick gold ring 
was inserted in the flange hole. Since the muon lifetime 
in gold is much shorter than in iron (tau ~ 0.073 ^s, 
Tpe ~ 0.2 ^s 1131), the time cut used during the analysis 
of the detected event substantially suppresses the back- 
ground arising from muon capture by the target body. 

Charged muon-capture products were detected by 
three silicon telescopes located directly in front of the 
kapton windows but still within the cooled vacuum envi- 
ronment (Si(7P, Sifl/, and Sioo in Fig.^. Each telescope 
consisted of two Si detectors: a 360 /xm thick dE/dx de- 
tector followed by a 4 mm thick E detector. The silicon 
detector preampliflers and amplifiers were RAL 108-A 
and 109, respectively |54j . Low-energy x rays from the 



muon cascade were detected by a 0.17 cm^ germanium 
detector (Ges in Fig. positioned outside the vacuum 
chamber, but separated only by several kapton windows 
from the reaction volume. Muon decay electrons were de- 
tected by four pairs of plastic scintillator counters (E^s, 
Eyp, Efl/, E/jo in Fig.^ placed around the target. 

The gas purity in the target was monitored by 75 cm'^ 
and 122 cm'^ germanium detectors [Ggm and Gep), 
which were sensitive to x rays between 100 keV and 
8 MeV. GeM and Ges were also used to monitor 
"harder" x rays, providing information about muon stops 
in the target walls. The NE213 detector was used to de- 
tect 2.5 MeV neutrons from dd fusion. 

The detector electronics triggering system was simi- 
lar to that used in experiments performed at TRIUMF 
(Vancouver, Canada) and details are given in Ref. fs^. 
The system measured events muon by muon, opening 
an 8 ^s gate for each received Tl pulse. At the end of 
the event gate, the individual detector electronics were 
checked and if any one detector triggered, all detectors 
were read and the data stored. If a second Tl signal ar- 
rived during the event gate, we assumed it was a second 
muon and discarded the event as pile-up. Great care 
was taken with the Tl threshold such that no muons 
would be missed, although this increased the rate of event 
gates started by electrons. Those events were rejected in 
software based on a lower-limit energy cut from the Tl 
scintillator. The pile-up rejection system was much im- 
proved over the TRIUMF version and reduced the detec- 
tion dead-time for multiple muons from w 50 ns down to 
3 ns. Thus we had only a 54 x 10^® chance per event to 
have two muons enter the target simultaneously without 
our awareness, although again an upper-limit cut on the 
Tl energy reduced the number of these events accepted 
in the analysis. 



B. Experimental conditions 



TABLE I: Experimental conditions. The last column, A^^, 
represents the number of muons stopped either in pure ^He 
or in the D2 + ''He mixture. 



Run 


Target 


Temp. 


Pressure 




C3Ho 






[K] 


[atm] 


[LHD] 


[%] [lO-^] 








6.92 


0.0363(7) 










6.85 


0.0359(7) 




I 


^He 


33 






100 1555.5 








6.78 


0.0355(7) 










6.43 


0.0337(7) 




II 


D2 + ^He 


32.8 


5.05 


0.0585(12) 


4.96(10) 4215.6 


III 


Da + ^He 


34.5 


12.04 


0.168(12) 


4.96(10) 2615.4 



The experiment was performed using three different 
gas conditions which are presented in Table The first 
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measurement, Run I, was performed with a pure ^He gas 
at different pressures. The second and third measure- 
ments used a D2 + ''He mixture at two different pres- 
sures. Run II was performed at 5 atm, whereas Run III 
took place at a pressure more than twice larger, namely 
12 atm, where it was necessary to raise the temperature 
to avoid liquefying the mixture. The density </> is given 
relative to the standard liquid hydrogen atomic density 
(LHD), iVo = 4.25 x 10^^ cm-^. As seen from the last 
column of Tabled Run II was by far the longest run be- 
cause its original purpose was to measure the fusion rate 
in the d/i^He molecule and the muon transfer rate from 
d/i atoms to '^He nuclei Isd]. 



A^ap, AJ?ap, and A*ap are the ^He total muon capture rates 
when producing a proton, Eq. a deuteron, Eq. 
and a triton, Eq. ||2J), respectively. An analogous equa- 
tion like Eq. Q should also be written for the production 
of deuterons. To avoid complication, we only write equa- 
tions for the protons using the p index. 

Thus the proton yield produced in the reaction 
during a time interval AT — [ti;t2] for the full energy 
range AE'™'*'^ is 



]\fHc 



'cap 



ft, 



(8) 



with the time factor ft given as 



III. MEASUREMENT METHOD 



ft = e-^Ho-«l . „ g-AHc-<5t^ 



(9) 



This section describes the method used to measure the 
differential muon capture rates by ^He nuclei with the 
production of protons and deuterons, as given in Eqs. (QJ 
and (O. Essentially, it is a simultaneous analysis of the 
time and energy spectra of events detected by the Si{dE— 
E) counters when muons stop in the gas target. 

The first step is to obtain time and energy spectra 
from the three Si{dE ~ E) detectors for each run. As a 
function of time, we then create two dimensional energy 
spectra {dE — E) to suppress essentially the accidental 
coincidence background and to separate precisely the two 
regions corresponding to the protons and deuterons. 

The second step is to simulate via Monte Carlo (MC) 
the time and energy distribution of the events detected by 
the 'B\{dE — E) detectors. The simulations are performed 
as a function of different proton and deuteron energy 
distributions. 

The final step is a comparison between the experimen- 
tal results and the MC simulation. The first comparison 
is done using the least-squares analysis between MC and 
data, and is described in Sec. lIV Al The second compar- 
ison requires first to transform the experimental spectra 
such that one obtains the initial energy distribution using 
Bayes theorem. This analysis is given in Sec. IIVBI 

The number of protons with a full kinetically allowed 
energy range AE™^^ — [0; E™^^] per unit of time is, for 
the case of pure ^He 



dNJAE'^'^'^^t) „ . , 

dt ^ P 



(5) 



where A^^° is the number of muons stopping in "^He and 
A^jjp is the muon capture rate in '^He when producing a 
proton. We use the rate Aho as the sum 



Ahc — An + a, 



He 
'cap 



(6) 



where Aq is the free muon decay rate (Aq — 0.4552 x 
10^ s~^), and Aj^^'p is the total muon capture rate in ^He, 
given by 



where 6t — t2 — ti (here and later in the paper we denote 
by Ax = [a;i;a;2] the interval of the quantity x and by 
5x = X2 — xi the interval width). 

We are now interested to know the proton yield for a 
certain energy range AEp — [Ep'.Ep + SEp] (the proton 



energy lies between Ep and Ep 
Np{AEp,AT) is then 

^He 



SEp). Such a yield. 



Np{AEp,AT) ^ -^-ft / 

^He J 



He Agap(Ai^p) .^^ (10) 
Ahc 



if one defines 



ALp(Ai?p) = j 



dXL 



By using Eq. (|10|) . one can write the capture rate as 
function of the energy range as: 



Xl^JAEp) = 



Np{AEp,AT) ■ Ahc 

N^'' ■ ft 
NpjAEp, AT^^^) ■ Xue 

ATHe 



(12) 



where AT™'^'^ = [0;oo). Therefore the differential cap- 
ture rate averaged over the proton energy range becomes 



dAg,p(gp) 
dE„ 



Np{AEp,AT) Ahc 1 
SEp ft 



(13) 



The number of muon stops in helium A'^'^ is found by 
measuring the yield and time distribution of muon decay 
electrons stopped in the target (gas and wall). The total 
number of muon stops is given by 



A, 



He 



\p I \d I \t 
^cap * ^cap ' ^cap " 



(7) 



/-wall 



(14) 
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The muon decay electron time spectra can be reproduced 
by a sum of exponential functions due to the muon stop- 
ping in aluminum and gold (target walls) as well as in 
the gas. 



dt 



where ^ai, ^Au, and Ahc are the normalization ampli- 
tudes and 



Al 



Aai — Qai • Ao -f- a 
Aau = Qau • Ao + A^jjp , 
Aho = Aq + Xf°p , 



(16) 



are the muon disappearance rates in the different ele- 
ments (the rates are the inverse of muon lifetimes in 
the target wall materials). In reality, Eq. (|15|l is an 
approximation of a more complex equation, which can 



be found in Ref. 



The nuclear capture rates in 



aluminum and gold, X^^^ = 0.7054(13) x lO*" s'^ and 
A^4, = 13.07(28) X 10^ s-\ are taken from Qai and 
Qau are the Huff factors, which take into account that 
muons are bound in the Is state of the respective nuclei 
when they decay. This factor is negligible for helium but 
necessary for aluminum Qai = 0.993 and important for 
gold Qau = 0.850 j53]. The constant B characterizes the 
random coincidence background. 
By measuring the amplitude Ajjo 



He 



Ee • Ao 



(17) 



and knowing the electron detection efficiencies Se aver- 
aged over the energy distributions, one obtains the num- 
ber of muons stopping in helium Nj^°. The muon decay 
electron detection efficiency Eg is determined experimen- 
tally as the ratio 



(18) 



where7Va;_e is the number of x rays of the /i'^He Ka line, 
measured by the germanium detector (Geg), in coinci- 
dence with a muon decay electron. is the same num- 
ber of X rays of the /i'^He Ka line when no coincidence is 
required. 

By determining the quantities Np{AEp,AT) and 
Nd{AEd, AT) based on the analysis of the two- 
dimensional energy distributions {dE — E) , knowing 
A^a°p = 2216 (70) s-i as determined in Ref. (this value 
is in good agreement with the calculated value A|?j[p = 
2140 s~^ from ^3), we can obtain the muon capture rate 
for protons X^^^{AEp) and deuterons X'^^^{AEd) as well 
as both differential rates dX^^p/dEp and dX'^^^/ dEd as a 
function of the proton (deuteron) energy. 



IV. THE ANALYSIS OF EXPERIMENTAL 
DATA 

As already mentioned in Sec. Ill Bl the experiment was 
performed using two different gases, namely pure '^He as 
well as a mixture of D2 -|- ^He. When a muon is stopped 
in the gas mixture, different processes occur. A diagram 
of processes occurring in the D2 -I- '^He mixture (the most 
complex one) is displayed in Fig. |21 

In the run with pure '^He (Run I of Table P) the quan- 
tities A^jjp(Ai?p) and AJ?j^p(A£'d) for the protons and the 
deuterons are found according to Eq. ((T^ . In the runs 
with a D2 -I- ^He mixture (Runs II and III of Table P) the 
same rates arc found as follows. The number of protons 
per time unit 



dNp 
IT 



^D-(e 

qis • aD • Ad3Ho • • C3 



■% 



(19) 



He 



A 



fid 



Ahc 



with Np^^" the number of muon stops in the D2 + '^He 
mixture and A^ajjo • • cajj^ the experimental molecular 
d/i'^He formation rate, using the known value A(j3He = 
2.42(18)xl0«s-i [33. TherateAHe = 0.457 xlO^s-Ms 
given by Eq. ^ using the known total capture rate ^3 ■ 
4> and cajje are the target density and helium concentra- 
tion given in Table P The experimental disappearance 
rate A^d for the d/j, atom in the ground state is given as 



A^d — Ao + Ad^d • ■ Cd ■ + Ad3 



He 



C3Hc 



(20) 



using the d/xd molecular formation rate Ad^d = 0.05 x 



|57l l58l l59l |60j and the effective muon sticking 



coefficient to the ^He nucleus resulting form the nuclear 
fusion reaction in the d/xd molecule d/zd /x^He -I- n, 
ujd — 0.07 [ssL l6l|. The deuterium concentration Cd = 
1 — C3He is obtained from Table D In reality, Eq. ^ 
is an approximation of a more complex equation, which 
can be found in Ref. |^ . 

The total probability ajje for /i'^He formation is 



aHc 



(21) 



where a^^ is the muon capture probability by '^He and 
ajjj, is the probability of muon transfer from an excited 
state of the d/i atom to '^He. Explicitly, 



aHo 



^Hc 



aD = 



A ■ C3He 
1 + A- C3He ' 
(1 - qis)aD , 
1 

1 + A- C3He ' 



(22) 



where A is the ratio between the stopping powers of the 
deuterium and helium atoms, A = 1.7 (2) and a^ is 
the muon capture probability by a deuterium atom, qis 



6 



[d^i^He]+ + e 
2p(J J=l 



a + M" +p(14.7MeV) 
+ 7(16.4 MeV) 




^Hc + + n(2.45 MeV) 
t + +p(3.02 MeV) 



\ j=o 



X'^=l \-^=i \'^^ 



[M^He + d] 



[d/i^He]+ + e 
2pcr J=0 



J 



\./=0 >J=0 \J=0 

P ' Ag 5 A^ 



FIG. 2: Scheme of muon processes in the D2 + ^He mixture. Muon capture by '^He occurs via process ^1 (with ~30% yield). 
The essential part of the capture (approximatively 65%) occurs after the d/i^He formation (#2). A small amount of capture is 
occurring after dfid fusion (#3). Details about all processes and rates are found in Ref. 



is the probability that the excited {dfj,)* atom will reach 
the ground-state. The term qis • aD is the probability 
for a muon stopped in the D2 + '^He mixture to be cap- 
tured by a deuterium atom and reach the ground-state. 
The qis values for the Runs II and III are 0.80 and 0.72, 
respectively, according to Refs. [63,E1j0I- These val- 
ues are somewhat higher than a^c~t ^periment M 
(qis — 0.689 (27)) performed at an intermediate '^He con- 
centration (csfjo = 9.13%). Using Eqs. and (|22J) one 
can then write 



1 



aHo 



1 + A- C3Hc 



(1 - qis + ^ • C,3Hc) . (23) 



Thus the proton yield in the time interval AT — [i 1 ; ^2] , 
for the whole energy range AE™'^^ is given by 



D/Ho 



with the time factor ft given as 



Ah. 



'cap n 



(24) 



-AHc-i2 



(25) 



The number of protons following muon capture in the 
energy range AEp is then 



NpiAEp,AT)^Nj^^^--ff 
and the capture rate becomes 



A?ap(A£;p) 



^He 



ALp(Ai?p) 



iVp(A£;p,Ar)-AHc 



D/Hc 



(26) 



(27) 



Note that Eqs. l(77|) and ifT^ are similar for both the 
pure gas and the mixture. The difference lies in the time 
factor ft, given by Eqs. and 

The calculation of ft for the D2 -I- '^He mixture 
(Eq. H25() ) demands the previous knowledge of ano, A(j3Hej 
Ad/id I Ahcj and A^d- Even if most of those values are 
known from other experiments, this experiment allows 
us another independent determination of these quantities 
and hence a consistency check. The rate of A^d is found 
by analyzing the time distribution of either the proton, 
the deuteron, or the photon emitted after d/i'^He forma- 
tion. The time distribution can be fitted using Eq. l(TI?|l . 
For Run I with pure '^He the value of ft was determined 



by using A^^^p in Eq. 



As mentioned in Sec. IIIII we find the number of muon 
stops in the gas by fitting Eq. (|15|) to the muon decay 
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FIG. 3: Time distribution of muon decay electrons measured 
in Run I with pure '^He. Top-right picture shows details of 
early times. 



electron time distributions. Figure |31 represents such a 
fit of electron time spectra when all four detector pairs 
E(7P, Eflj, EuOj and Fj^e are added together. 

Figure 21 displays the energy spectra of the low-energy 
photons from /x^He atoms (Ka at 8.2 keV, K/3 at 9.6 keV, 
and K7 at 10.2 keV) measured with the germanium de- 
tector Ges with and without the delayed muon decay 
electron coincidence. The electron detection efficiency 
Se is determined using Eq. IjlSI) . The so obtained value 
still needs to be corrected for the difference in posi- 
tions between the germanium and the Si{dE — E) de- 
tectors with respect to the muon stop distribution along 
the incident muon beam. The final value for the to- 
tal muon decay electron detection efficiency of the four 
electron counters found from the analysis of Run II is 
Ee = 16.4 ±0.22% IHIi^. 

Since the background is mainly caused by muon stops 
in the target walls (Al, Au) followed by their nuclear cap- 
ture and the emission of charged products (with charac- 
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FIG. 4: Muonic x-ray spectra measured by the germanium 
detector in a mixture of D2 -I- "^He, without (open circles) and 
with (solid circles) muon decay electron coincidences. 



teristic times tai = 0.865 ^s and tau = 0.073 fj,s l53j). 
the background contribution will be determined in two 
steps. 

The first step is to remove the background contribution 
from muon stops in gold. Hence, we selected only events 
detected by the Si{dE — E) detectors for times t > Atau- 
The remaining events are due to muon stops in the gas, 
which have a time distribution following Eq. (O for pure 
■^He and Eq. H19() for the mixture D2 -I- '^He, and muon 
stops in aluminum. Therefore, the time distribution of 
the Si{dE — E) events in pure ^He is represented as: 



Jt 



Di ■ e" 



D2 



c 



with 



Di = 
D2 = 



Al 



vAl 



~A1 
P 



•ALp(Ai?p) 



-He 



(28) 



(29) 



^ . ^ represent the number of muons 



The terms Nj^"" and iV,^' 
stopping in helium and aluminum, respectively, ' and 
are the proton detection efficiencies after muon cap- 
ture in aluminum or helium averaged over the energy 
interval AEp, and C is the accidental coincidence back- 
ground. 
For the D2 - 

as 



^He mixture, Eq. H28|) has to be rewritten 



di 



= Di-e 



— AArt_ 



'Die 



-A„ 



-C (30) 



with 



D' 



7VD/Hc.AP^p(Ai?p).£f .(aHc + &) , 
</H^-A?,p(Ai?p).£H°.eD , (31) 



where the term nJ^^^'^ represents the number of muons 
stopping in the mixture. The constant D2 is replaced by 
the corresponding D'2 and D3. 

Figure |S1 displays the time distributions of Si{dE — E) 
events for the experiment with pure '^He. The time distri- 
butions are very well fitted by Eq. (|28|) . using the values 
Aai = 1.156 X lO^s"^ and Ahc = 0.457 x lO^s^^ from 
Refs. 

The second step is to remove the background arising 
from muon stops in aluminum. For this purpose, the 
time interval St — t2 — <i is divided into two subintervals 
6tA — — ti and Sts = t2 — t-^. Therefore the proton 
yields Np which correspond to the two new intervals 6tA 
and 6tB have the form 



dN, 



dt 



^dt 



Aai 
D2 



-Aai-*! 



'l-e 



-AAr<5*A ^ 



(32) 



+ 2^_AHc-tl(l_g-AHo-5tA)_^C'.J^^ 

A Ho 
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FIG. 5: Time distributions of Si(di5 — E) events for Run I: 
(a) protons, (b) deuterons. The histograms represent the ex- 
perimental data. The solid lines 1 and 2 are the exponential 
functions for Al and '^He, whereas 3 is the accidental back- 
ground coincidence. 



and 



dt 



Aai 

+ iZle"^«=*^(l-e-^«<='^'*«)-HC-fts. 

Aho 



(33) 



The total numbers of events and A'^f , given for the 
two time intervals SIa and Sts are given by the two- 
dimensional amplitude distributions {Aj-k)A and (yljfe)^, 
created for each (jk) cell, where j = 1, . . . ,^ and fc = 
1, . . . , m are the cell indexes on the dE (the energy losses 
in the thin Si detector) and the E axes (the deposited 
energy in the thick Si detector), respectively. 

The time intervals, 8tA and ^^b, are chosen such that 
the difference between the proton or deuteron yields mea- 
sured in the intervals 6tA and btB is independent of the 
aluminum muon capture contribution. This means that 
the first parts of Eqs. H32|) and H33|l are then equal, i.e.. 



Aai 



^1 



) 



Aai 



(34) 



If the initial t\ and final ti measurement times are 
given, the middle time becomes 



Aa: 



■ ln2. 



(35) 



The difference between and iV^ is the total num- 
ber of events in the resulting {^Ajk)B-A two-dimensional 
(dE — E) protons distribution. This distribution was ob- 
tained by subtracting channel by channel the content of 
the (jfc) cell for the two (^Ajk)A and (^Ajk)B distributions. 

The final number of protons N^''^'^^ for the pure '^He 
measurement is then 



final 



with 



C ■ [5ti 



g-AHo-tl(^ 
g-AHo-t3(l 



cap V 
Ahc 
5tA] 



(36) 



-Aho-i5*£ 



whereas, for the D2 -I- '^He mixture, it becomes 



rfinal 

p ' '-'P ' 'cap \ 

-I- C ■ {btB - ^tA\ 



(37) 



(38) 



with 



Ft 



(aHc + Cd) 



2g~'^Hc-i3 g — Ahc'^I g — AHc-i2 



A 



Ho 



(39) 



g — A,id-*i^]^ _ g— A^d-i5tA^ 



A 



Analyzing the data according to Eqs. H34(l and H35(l 
we obtained the intervals bdA — [ii;i3] = [0.51; 1.098] ^s 
and = ^z^t-^ = [1.098; 6.0] ^s. The corresponding 
capture events in aluminum amount to ~ 23% of the 
total events. As an example. Table HH show the number 
of events measured in Run II in both time intervals and 
both elements. 

Our subtraction method, while reducing the number 
of events in helium by a factor 2 (see Table , yields 
essentially background-free events. However, Eqs. ifll^ 



TABLE 11: Number of aluminum and helium events in Run II, 
as a function of the time intervals AtA and Ats. 



Particle 


Interval IStA 


Interval Ats 




Aluminum Helium 


Aluminum Helium 


Proton 


2700 3600 


2700 10 800 


Deuteron 


1150 1650 


1150 5800 
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FIG. 6: Two-dimensional energy distributions of the Si(d£' — 
E) detector events for the time interval 0.5 < t < 6.0 fis rela- 
tive to the muon stop, (a) represents Run I (pure ^He) and (b) 
Run II (D2 + ''He mixture). The distinction between protons 
and deuterons is clearly visible for both measurements. 



and (|38|l still contain some parameters that need to be 
determined, namely the energy interval AEp and the ac- 
cidental coincidence background described by the con- 
stant C. 

The energy intervals for detecting protons and 
deuterons by the Si{dE — E) detectors were chosen such 
that the real detection sensitivity is the same for any ini- 
tial energies. This allows us to remove any possible dis- 
tortion in our amplitude distribution, which would occur 
for too low or too high energies. The chosen limits are 
4 — 23 MeV for both protons and deuterons in the thick 
E detector. The thin dE detector has two different en- 



ergy intervals, namely 1 — 6 MeV for the protons and 
2-8 MeV for the deuterons. 

Figure El displays the two-dimensional [dE ~ E) distri- 
butions of events detected by the Si{dE — E) detectors 
in Run I with pure ^He and in Run II with the D2 -I- '^He 
mixture. The two distinct branches of events correspond- 
ing to the protons and the deuterons are clearly visible 
and lie inside our chosen energy intervals. Note that the 
shapes of the two-dimensional {dE — E) distributions ob- 
tained in the runs with pure '^He and with the D2 -|- '^He 
mixture coincide. This indicates that there are no ne- 
glected systematic errors, and that the algorithm used 
for the data analysis is correct. 

As to the accidental coincidence background described 
by the constant C, its contribution to Eqs. 13()|) and 138|) 
is small when compared to the muon stop contributions 
in Al and ^He, as can be seen in Fig. [S] The constant C 
was quantitatively determined in each Run by fitting the 
time distribution, as given in Fig. [S] including the time 
interval —0.4 /its < t < with respect to the muon stop. 
Details of such a fit shown in the muon decay electron 
time spectra are in Fig. |3| 

As mentioned in the introduction, we want to deter- 
mine different characteristics of the muon capture by '^He 
nuclei, namely: 

- the initial energy distributions of protons and 
deuterons {S{Ep), S{Ed)), 

- the muon capture rates as function of the energy 
for both the protons and deuterons {XP^p{AEp) and 
Aip(Ai?rf)), 

- their derivatives dX^^^/dEp and dXf^^/ dEd- 

For this purpose, following Eqs. ((T^ . (|77|l . and ((T^ . we 

need to determine the number of protons Np{AEp, AT) 
and deuterons Nd{AEd, AT) for each energy interval 
AEp and AEd- In the next two subsections we describe 
the two approaches to determine the respective number 
of protons and deuterons, based on the analysis of the 
two-dimensional {Ajk)B-A distributions as function of 
dE and E for each of the three Runs (I-HI) . 

A. Method I: Least— squares 

The principle of this method is to use Monte Carlo 
(MC) simulations to reproduce the experimental data 
and to minimize the free parameters which are required 
by such a simulation. The simulation conditions and pa- 
rameters will be given below. The energy spectra of the 
protons and deuterons produced by reactions and (0) 
are divided into i subintervals of 1 MeV fixed widths. 
Since the theoretical maximum energies are » 53 MeV 
for the protons and ^ 33 MeV for the deuterons, the 
numbers of subintervals are 53 and 33, respectively. 

Using the experimental muon stop distribution in our 
target, we simulate the probability P^'~^{Ajk/ E^) that 



10 



a proton (analogously a deuteron) produced with an en- 
ergy Ep (in the i-th interval AE^) will be detected by 
the Si{dE — E) detectors in the (jfc) cell of the two- 
dimensional distribution Ajk- This probability is 



MC 



ino)i 



(40) 



where (njfe)f"^ is the number of simulated events de- 
tected in the (jfc) cell when the number of protons, which 
were created with an initial energy E^ from the interval 
AEp, is {No)i. The {jk) cell size is chosen arbitrarily 
and mainly depends on the statistics of events (?^jfe)f"^ 
belonging to a particular cell (jk). 

Then the MC simulated "pseudo-experimental" (i.e., 
normalized to the experimental counts A^^"*' = 
Np{AEp, AT)) event numbers {Njk)^'^ for each (jk) cell, 
become 



i 

p 

S{Ep)dEp , (41) 



where S{Ep) is the initial proton energy distribution nor- 
malized to unity in the full energy interval AEp. 

In our energy intervals 10 < _Ep < 49 MeV and 
13 < i?d < 31 MeV both the proton and deuteron energy 
distributions, S{Ep) and S'(i?d), obtained via the impulse 
approximation model and the realistic wave functions for 
the ^He nucleus ground state [4^^|4l,|53 can be weU 
described by the expression 



SiEp) 



ApE 



-a„E„ 



where the amplitude Ap and the fall-off yield ap are the 
variable parameters. Thus Eq. 141|) can be rewritten as 

(Njk)''^ = iV«"-' J2 P^'^^iAjk/E;) • S{El,) , (43) 

i 

where 

s{e;) = s{e;) . (44) 



We created [Njk) for different values of the ampli- 
tude Ap and the fall-off yield ap and used the min- 
imization procedure between the MC and experimental 
events, 



EE 

]=i k=i 



MCI 



(45) 



to obtain the best values for the parameters Ap and ap 
which describe the initial energy distribution of protons 
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FIG. 7: E axis projections of the two-dimensional {dE — 
E) distributions for the protons (a) and the deuterons (b) 
obtained in Run II with the D2 -l-'^He mixture. The histogram 
shows the experimental data, whereas the black triangles are 
MC events from Method I. 



S{Ep). {Njk)^^^ is the number of measured events be- 
longing to the [jk) cell, obtained for pure ^He and the 
D2 + '^He mixture, respectively. The A'^"^' values of 
(42) Eqs. and represent the sum of the {NjkY'^^ over 



i and k. 

A second and parallel minimization is done when pro- 
jecting the experimental and MC events onto the two 
energy axes j and k. When projecting onto the E axis, 
we have the experimental data as 



(46) 



and the MC events as 



{Nk) 



MC 



= ^p"'' E E P'^'^iA^k/Ei) ■ s{e;) 

Therefore becomes 



(47) 



^ [{Nk) 



cxp 



(Nk 



iMCl 



j2 



(48) 



11 



Similar equations can be written for the second axis j 
when we project the events onto the dE axis. 

Figure displays the least-squares comparison of the 
E axis projection of the two-dimensional experimental 
and the MC simulated distributions for the protons and 
the deuterons of Run II. As seen, the MC distributions 
correspond very well to the experimental proton and 
deuteron energy distributions, thus strongly supporting 
our analysis method I. 

The amplitude and fall-off yield results from the three 
experimental Runs (I-III) are 

Ap = (0.832 ± 0.043) MeV"^ 

ap = (-0.163 ± 0.002) MeV~^ (49) 

for the protons and 

Ad = (5.59 ± 1.39) MeV"^ 

ad = (-0.243 ± 0.012) MeV"^ (50) 

for the deuterons. 

The capture rates XP^p{AEp) are obtained after using 
Eq. to calculate the proton yield Np{AEp, AT) and 
then applying Eqs. (|12|) and H27I) . The differential cap- 
ture rates dX^^^/dEp also follow from the proton yield 
and Eq. H13|) : they are given in Figs. 1121 and 1131 for the 
protons and deuterons, respectively. 
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FIG. 8: Experimental energy distributions (open triangles) 
of protons (a) and deuterons (b) found by the method I and 
averaged over Runs (I-III) in comparison with the theoretical 
model |4^ (solid line). 



TABLE III: Mean proton energy distribution normalized to 
unity in the energy range 10 < Ep < 49 MeV, from methods 
I and II. 





< S{Ep) > [McV- 


-^1 


[McV] 


Method I 


Method II 


10.5 


0.150 (13) 


0.1570 (83) 


11.5 


0.127 (12) 


0.1309 (48) 


12.5 


0.108 (10) 


0.1077 (35) 


13.5 


0.0922 (88) 


0.0958 (29) 


14.5 


0.0784 (77) 


0.0765 (24) 


15.5 


0.0667 (67) 


0.0644 (22) 


16.5 


0.0568 (59) 


0.0525 (22) 


17.5 


0.0483 (51) 


0.0485 (20) 


18.5 


0.0411 (45) 


0.0392 (18) 


19.5 


0.0349 (39) 


0.0313 (17) 


20.5 


0.0297 (34) 


0.0284 (16) 


21.5 


0.0252 (30) 


0.0251 (14) 


22.5 


0.0215 (26) 


0.0208 (14) 


23.5 


0.0183 (22) 


0.0184 (13) 


24.5 


0.0155 (20) 


0.0162 (11) 


25.5 


0.0132 (17) 


0.0150 (12) 


26.5 


0.0112 (15) 


0.01135 (28) 


27.5 


0.0095 (13) 


0.00934 (19) 


28.5 


0.0081 (11) 


0.00793 (16) 


29.5 


0.00688 (98) 


0.00679 (14) 


30.5 


0.00585 (85) 


0.00585 (22) 


31.5 


0.00497 (74) 


0.00528 (27) 


32.5 


0.00422 (64) 


0.00440 (21) 


33.5 


0.00359 (56) 


0.00379 (16) 


34.5 


0.00305 (49) 


0.00316 (13) 


35.5 


0.00259 (42) 


0.00255 (11) 


36.5 


0.00220 (37) 


0.00210 (9) 


37.5 


0.00187 (32) 


0.00177 (7) 


38.5 


0.00158 (28) 


0.00142 (6) 


39.5 


0.00135 (24) 


0.00119 (5) 


40.5 


0.00114 (21) 


0.00105 (4) 


41.5 


0.00097 (18) 


0.00092 (4) 


42.5 


0.00082 (16) 


0.00079 (3) 


43.5 


0.00070 (13) 


0.00067 (3) 


44.5 


0.00059 (12) 


0.00057 (2) 


45.5 


0.00050 (10) 


0.00048 (2) 


46.5 


0.00043 (9) 


0.00041 (2) 


47.5 


0.00036 (8) 


0.00034 (1) 


48.5 


0.00031 (7) 


0.00029 (1) 



The average energy distributions < S{Ep > and < 
S{Ed) > from Runs (I-III) normalized to unity for the 
energy intervals 10 < < 49 MeV and 13 < < 
31 MeV are given in Table IIIII for the protons and in 
Table IIVI for the deuterons. Figure |H1 displays the en- 
ergy distributions < S{Ep) > and < S{Ed) > averaged 
over Runs (I-III) in comparison with the model distri- 
butions obtained when treating the muon capture in the 
simple plane-wave impulse approximation |45| and in the 
impulse approximation with the realistic Bonn B poten- 
tial (621 of NN interaction in the final state . 

Experimental and theoretical results agree quite well 
within the statistical errors for the energy ranges 10 < 
Ep < 40 MeV and 13 < Ed < 24 MeV, respectively. 
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TABLE IV: Mean deuteron energy distribution normalized to 
unity in the energy range 13 < Ed < 31 MeV, from methods 



I and II. 


Sd 


< S(Ea} 


> [McV^^] 


[McV] 


Method I 


Method II 


13.5 


0.210 (44) 


0.216 (11) 


14.5 


0.167 (36) 


0.1690 (65) 


15.5 


0.133 (29) 


0.1281 (47) 


16.5 


0.106 (24) 


0.1043 (41) 


17.5 


0.084 (19) 


0.0842 (35) 


18.5 


0.067 (16) 


0.0674 (29) 


19.5 


0.053 (13) 


0.0521 (25) 


20.5 


0.042 (10) 


0.0426 (23) 


21.5 


0.0328 (84) 


0.0345 (23) 


22.5 


0.0258 (68) 


0.0251 (21) 


23.5 


0.0202 (55) 


0.0181 (18) 


24.5 


0.0158 (44) 


0.0132 (18) 


25.5 


0.0124 (35) 


0.0124 (17) 


26.5 


0.0096 (28) 


0.0101 (16) 


27.5 


0.0075 (23) 


0.0071 (16) 


28.5 


0.0058 (18) 


0.0058 (19) 


29.5 


0.0045 (14) 


0.0052 (18) 


30.5 


0.0034 (12) 


0.0044 (14) 



For proton energies Ep > 40 MeV and deuteron ener- 
gies Ed > 24 MeV a discrepancy exceeding the tolerable 
range determined by the statistical errors is observed. 
The cause of the discrepancy is not clear yet. It may 
be due to the necessity of taking into account exchange 
current contributions in the interaction and nucleon pair 
correlations in muon capture by the '^He nucleus. 



B. Method II: Bayes theorem 

In this approach we use the Bayes theorem |6^ [t^, 
0,[72 to determine the initial energy distribution, S{E), 
of the protons and the deuterons produced by muon cap- 
ture in ■^He. For this purpose, we apply inverse trans- 
formations from the detected two dimensional {dE — E) 
amplitude distributions. 

The relation between the probability P{Ajk/ E^) that 
a proton produced with an initial energy E^ (in the i- 
th interval of 1 MeV width in our case) will be detected 
by the Si(d£' — E) telescopes and the inverse probability 
P{Ep/Ajk) (probability that a proton detected in the 
{jk) cell comes from the AE^ subinterval) is 



PiE;/A,,) - 



s{e;) ■ P{A,k/E;) 
T.~S{El,)-P{A,k/El) 



(51) 



The probability P{Ajk/Ep) is given by the MC simulated 
probability P"''^ {A^k/ E^) defined in Eq. gDJ. 

In the first step of the analysis we start from the initial 
energy distribution So{Ep) = S{Ep) given by with 
an arbitrary set of parameters. When using the probabil- 
ity given by Eq. (|51|l and the experimental data of each 
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FIG. 9; E axis projections of the two-dimensional {dE — E) 
distributions for the protons (a) and deuterons (b) obtained in 
Run II. The histogram shows the experimental data, whereas 
the black triangles are MC events from Method II. 



(jk) cell, we obtain a set of i relations 



Np{AEp,AT)-S{El)^ 



E E PiE;/A,k) 
j=i k=i 



PiA/E;) 

(52) 

where Np{AEp, AT) corresponds to Eqs. ifTU)! or if^ . 
and P{A/Ep) is the probability that a proton of initial 
energy Ep is detected anywhere in the proton branch of 
the two-dimensional distribution Ajk- This probability 
can be written as 



PiA/E;)^Y.J2P 

J=l fc= 



MC 



{A.k/E;) 



(53) 



We then compare Np{AEp, AT) and the experimental 
counts A'^'^'^P = J2J2i^3k)°^''^ for each i-th interval via 
a analysis and obtain a proton energy distribution 
S{E^) from Eqs. ((S3l and As long as the is 

not satisfactory, we re-use the last S{Ep) as the starting 
values in Eq. H51() in the next iteration. 

In addition, the initial energy distributions of the pro- 
tons and deuterons can also be derived by analyzing 
the projections of the two-dimensional distribution (Aj^) 
onto the dE axis (Aj) and the E axis (Ak). The equa- 
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FIG. 10: Experimental energy distributions (open triangles) 
of protons (a) and deuterons (b) found by the method II and 
averaged over Runs (I-III) in comparison with the theoretical 
model (solid line). 



tions for the dE axis are 



P{E\,IA,) = 





fe=i 



(54) 



and 



P{AIEI 



(55) 



Similar equations can be written for the E axis. Using 
the above equations, we obtain simulated values for the 
proton and deuteron yields as measured by the Si{dE—E) 
detectors. Figure El shows the projections of the experi- 
mental and simulated {dE — E) distributions for protons 
and deuterons onto the E axis. 

The mean proton < S{Ep) > and deuteron < S{Ed) > 
energy distributions from Runs (I-III) are given in Ta- 
bles IHII and CxI The mean values are also displayed in 
Fig. ^1 It is important to note that the distribution 
S{Ep) does practically not depend on the form of the en- 
ergy distribution So{Ep) which is chosen for the first iter- 
ation. Variation errors in the determination of S{Ep) fall 
within the statistical errors of {Nj^Y^^^ ■ Since Eqs. H52|) 
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FIG. 11; Comparison of the proton (a) and deuteron (b) en- 
ergy distributions found by methods I (black triangles) and II 
(open circles) averaged over Runs (I-III). For sake of visibility, 
we plotted both methods results alternatively. 



and (|55|l (as well as the other projection) have an iden- 
tical solution, their comparison makes it possible to con- 
clude, with an accuracy determined by the statistics of 
the detected events, that there are no systematic errors 
in the analysis of experimental data. 

A comparison between the experimental energy dis- 
tributions given in Fig. 1101 and the energy distribution 
calculated by the impulse approximation reveals some 
discrepancies of the same character as in the method I 
analysis, as long as the interactions between the reaction 
products 1^ and (0) are considered and a realistic Bonn 
B 48] nucleon-nucleon potential is employed. 

The capture rates \^^^^{/S.Ep) as well as the differential 
capture rates dX^^^/dEp which are found from Eqs. H52|) 
and ifS^ using Eqs. ifT^ . ifT^ . and (|77|) are given in 
Fig- El for the protons and in Fig. El for the deuterons. 



V. CONCLUSIONS 

The proton and deuteron energy distributions found 
by methods 1 and II largely coincide within the mea- 
surement errors, which points to the compatibility of the 
different approaches and to the absence of any systematic 
errors which may have been neglected in the analysis of 
the experimental data (see Fig. lll|) . However, the errors 
on S{Ep) and S{Ed) found by both methods are differ- 
ent. The analysis using method 11 gives a more precise 
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TABLE V: Relations d\^^p{Ep) /dEp found by methods I and 
II and averaged over Runs (I— III). The proton energies Ep 
correspond to the middle of the respective 1 MeV intervals. 



Ep 
[MeV] 


Mcthod I 


Method II 


10.5 


5.49 (59) 


5.77 (47) 


11.5 


4.67 (51) 


4.81 (35) 


12.5 


3.98 (44) 


3.95 (28) 


13.5 


3.38 (38) 


3.52 (25) 


14.5 


2.88 (33) 


2.81 (20) 


15.5 


2.45 (29) 


2.37 (17) 


16.5 


2.08 (25) 


1.93 (15) 


17.5 


1.77 (22) 


1.78 (13) 


18.5 


1.51 (19) 


1.44 (11) 


19.5 


1.28 (16) 


1.151 (95) 


20.5 


1.09 (14) 


1.041 (88) 


21.5 


0.93 (12) 


0.920 (77) 


22.5 


0.79 (11) 


0.763 (71) 


23.5 


0.671 (92) 


0.675 (64) 


24.5 


0.570 (80) 


0.595 (56) 


25.5 


0.485 (69) 


0.549 (55) 


26.5 


0.412 (60) 


0.417 (28) 


27.5 


0.350 (52) 


0.343 (23) 


28.5 


0.298 (45) 


0.291 (19) 


29.5 


0.253 (39) 


0.249 (17) 


30.5 


0.215 (34) 


0.215 (16) 


31.5 


0.183 (29) 


0.194 (16) 


32.5 


0.155 (26) 


0.162 (13) 


33.5 


0.132 (22) 


0.139 (11) 


34.5 


0.112 (19) 


0.116 (9) 


35.5 


0.095 (17) 


0.094 (7) 


36.5 


0.081 (14) 


0.077 (6) 


37.5 


0.069 (12) 


0.065 (5) 


38.5 


0.058 (11) 


0.052 (4) 


39.5 


0.050 (9) 


0.044 (3) 


40.5 


0.042 (8) 


0.038 (3) 


41.5 


0.036 (7) 


0.034 (3) 


42.5 


0.030 (6) 


0.029 (2) 


43.5 


0.026 (5) 


0.024 (2) 


44.5 


0.022 (4) 


0.021 (2) 


45.5 


0.019 (4) 


0.018 (1) 


46.5 


0.016 (3) 


0.015 (1) 


47.5 


0.013 (3) 


0.013 (1) 


48.5 


0.011 (2) 


0.011 (1) 



information about the proton and deuteron energy dis- 
tributions than method I. In method I, we compare using 
the numbers of detected events from a (jfc) cell with sim- 
ilar MC simulated data. Such numbers are the sums of 
the contributions from all z-th proton energy subintervals 
AEp. In method II, we have much deeper relations be- 
cause the comparisons are performed via Eq. (|52() for each 
i-th subinterval separately and all comparisons should be 
simultaneously satisfactory. 

Similar remarks hold for the differential capture rates 
dXP^jEp)/dEp and \P^p{AEp), as seen in Fig. [El Ta- 
bles El and |^ list the values of < dXP {Ep)/dEp > 



TABLE VI: Relations dAcap(-Bd) /dEa found by methods I and 
II and averaged over Runs (I-III). The deuteron energies Ed 



correspond to the middle of the respective 1 MeV intervals. 


E^ 
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23.5 


0.43 (12) 


0.397 (47) 


24.5 


0.340 (95) 


0.289 (44) 


25.5 


0.266 (76) 


0.272 (41) 


26.5 


0.207 (61) 


0.221 (37) 


27.5 


0.161 (49) 


0.156 (35) 


28.5 


0.124 (39) 


0.127 (42) 


29.5 


0.096 (31) 


0.114 (41) 


30.5 


0.074 (25) 


0.095 (31) 



and < dX^^p{Ed) / dEd > found from the analysis of the 
Runs (I-III) with pure '^He and D2 + '^He mixtures data 
by methods I and II. 

The addition of the differential rates in Tables 
and IVII yields the muon capture rates by the ^He nu- 
cleus followed by proton and deuteron production in the 
final state in the energy intervals 10 < Ep < 49 MeV and 
13 < -Ed < 31 MeV, respectively (see Table Ivin . 

Looking more closely at Figs.Eland^l our results and 
their comparison with the experimental data [3^l40l | and 
the calculations [isl liM l50l | indicate the following results 
for the protons and deuterons. Experimental (obtained 
by methods I and II) and calculated differential rates 
d\P^^{Ep) / dEp and XP^p{AEp) for muon capture by the 
■^He nucleus followed by proton production in the energy 
range 10 < Ep < 40 MeV show quite good agreement 
both in form and in magnitude. The calculations were 
carried out in the simple plane-wave impulse approxima- 



TABLE VII: Muon capture rates by ^He nucleus (in s"^) fol- 
lowed by the proton and deuteron production following meth- 
ods I and II. 



Method 


I 


II 


A?ap(10<i5p<49MeV) 


36.7 ± 1.2 


36.8 ±0.8 


Alp(13<i?p<31MeV) 


21.3 ± 1.6 


21.9 ±0.6 
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FIG. 12: Differential rates d\^^p{Ep)/dEp (open circles) 
found by methods I (a) and II (b) averaged over Runs (I- 
III). Black triangles are the results of Refs. pffll. l4(il|: the solid 
line corresponds to the model j45|. 
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FIG. 13: Differential rates d\'^^p{Ed)/dEd (black triangles) 
found by methods I (a) and II (b) and averaged over Runs (I- 
III). Black boxes are the results of Refs. |39l |40| : the solid 
line corresponds to the model |45|: the dotted line is based on 
calculations from Ref. |4^ . 



tion (PWIA) with allowance made for final-state interac- 
tion of reaction ^ and ^ products. However, there is 
a difference between the results of the present paper and 
the calculations for proton energies Ep > 40 MeV. 

The measured dependence dX^^p{Eii)/dEci found by 
using methods I and II is quite well described by the 
theoretical PWIA dependence 01 in the deuteron en- 
ergy ranges 13 < Ed < 20 MeV (Method I) and 13 < 
Ed < 17 MeV (Method II), respectively. For deuteron 
energies Ed > 20 MeV there is a noticeable discrepancy 
between experiment and theory The measured val- 
ues of dX'^^p{Ed)/dEd and the PWIA calculations ^ 
with the refined realistic NN interaction potential (Bonn 
B) appreciably disagree over the entire deuteron energy 
range. 

Next, we can estimate the total capture rate (full en- 
ergy range [0; oo)) using a simple extrapolation of our 
data at low energies and a one-exponential weighted fit 
of the differential capture rate in the full energy range. 
Using the function 

^^^(EA = H ■ e-^ ''^ , (56) 
dEp 

where H and G are free parameters, we obtain the total 
capture rate for the proton as their ratio 



Results for protons and deuterons, using both meth- 
ods, are given in Table Ivml The summed rate A^^^p -I- 
•^cap (which corresponds to Eq. O without the tri- 
ton contribution) is also compared to other experimen- 
tal 41, 42, 42J and theoretical ^,^^J values. Agree- 
ment between our results and previous ones is excellent. 

An experimental determination of muon capture on 
■^He nuclei makes a study of electromagnetic and weak 
interactions of elementary particles with 3N systems pos- 
sible without introducing uncertainties due to inadequate 



TABLE VIII: Total muon capture rate for reactions Q 
and The results of this work is an estimation from both 
methods I (least-square) and II (Bayes theorem). 



Method 


\p 


'^cap 


\P 1 \d 
'^cap ' ''cap 










This work 








Least-square 


187 ± 11 


491 ± 125 


678 ± 126 


Bayes 


190 ±7 


497 ± 57 


687 ± 60 


Zaimidoroga 






660 ± 160 


Auerbach [42j 






665 li™ 


Maev [43] 






720 ± 70 


Yano [44] 






670 


Philips [45] 


209 


414 


623 


Congleton [46] 






650 
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approximations of 3N states in the analysis. According to 
the theory, meson exchange currents must also be taken 
into account in future analysis of experimental data. As 
compared with Refs. 39, 40], this experiment yields for 
the first time information on the "softer" region of proton 
and deuteron energy spectra, which is more sensitive to 
the theoretical models describing the final-state nucleon- 
nucleon interactions. 

Finally it should be mentioned that by increasing the 
efficiencies of the proton and deuteron detection systems 
and their functional capabilities, by decreasing the lower 
and increasing the upper thresholds in the Si{dE — E) 
telescopes, the above method will provide precise infor- 
mation on the characteristics of muon capture by bound 
few-nucleon systems. It then becomes possible to ver- 
ify various theoretical models of muon capture by helium 



nuclei and to clarify the nature of discrepancies between 
the results of the present paper and the experimental 
data lailiS- 
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